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Adhesive contact between a thin elastic sheet and a substrate in a liquid environment arises 
in a range of biological, physical and technological applications. By considering the dynamics of 
this process that naturally couples fluid flow, long wavelength elastic deformations and microscopic 
adhesion, and solving the resulting partial differential equation numerically, we uncover the short- 
time dynamics of the onset of adhesion and the long-time dynamics of a steady propagating adhesion 
front. Simple scaling laws corroborate our results for characteristic waiting-time for adhesive contact, 
as well as the speed of the adhesion front. A similarity analysis of the governing partial differential 
equation further allows us to determine the shape of a fluid filled bump ahead of the adhesion 
zone. Finally, our analysis yields the boundary conditions for the apparent elastohydrodynamic 
contact line, generalizing the well known conditions for static elastic contact while highlighting how 
microscale physics regularizes the dynamics of contact. 


INTRODUCTION 

Adhesion of thin elastic films to a substrate occurs in a range of problems in physical chemistry [T] , in biology [2] , 
and in engineering [SHin]. Experimental and theoretical studies of adhesion focus on the microscale physical chemistry 
of bond formation mini, and the mesoscale mechanics of elastic deformation of the film, and squeeze flow in the 
gap IT 2 H 2 T]. In the neighborhood of the contact line itself, all these effects are potentially important - short range 
adhesive interactions affect the squeeze flow, which can lead to a high pressure that deforms the film and affects the 
dynamics of adhesion in return. However, theoretical models of soft adhesion usually focus only on the elastic and 
adhesive aspects of the process, typically neglecting the hydrodynamics [22II24] . An important exception is [8] where 
energetic arguments were used to determine the steady bonding speed of silicon wafers coming into contact and the 
self-similar shape of the contact zone, but without a complete consideration of the microscopic nature of adhesion or 
the transient dynamics of contact. 

Here we consider the complete transient elastohydrodynamics of soft adherent sheets with the goal of characteriz¬ 
ing the onset of adhesion, the generic properties of the moving contact zone once it reaches a steady state. We also 
determine the appropriate boundary conditions at the effective contact line (Fig. 1), thus generalizing the bound¬ 
ary conditions for a static elastic contact line, first studied theoretically and experimentally nearly a century ago 
[28] . These questions are analogous to those raised a long time ago for liquid-vapor-solid contact lines, which have 
an accompanying rich literature [25ti27] . Just as in liquid-vapor-solid contact lines, elastohydrodynamic contact is 
associated with a near-singular mathematical behavior. 



FIG. 1: Sketch showing variables. An elastic sheet with a thickness b, which is separated from a solid substrate by a thin 
viscous film of height h = h{x,t). The sheet is adhered to the substrate at the left side at a height h — a and has at the right 
side a far-field height /loo- The adhesion pressure brings the sheet towards the substrate leading to the formation of an 
adhesion zone R{t) that propagates towards the right. 
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FIG. 2: Elastohydrodynamics of adhesion. A thin and long elastic sheet that is initially separated from a solid substrate by 
a viscous him is brought in contact at the left edge by a microscopic adhesion potential. The shape of the sheet follows from 
the evolution equation Q, starting with an initial condition given as h{x,i = 0) = 1 + 0.1 x (1 — tanh(25x)^) that trigger 
adhesion. The results are presented in non-dimensional form for hve different points in time t = [tw = 47,100,200,300,400], 
with B = 5 X 10“^, e = 10“^ and a = 0.12. After a short transient phase associated with the formation of a contact zone 
given by the waiting-time Tw = — e“^( 3 B) 2 , the adhesion zone propagates in time displacing the huid into a bump. 

The position of the apparent elastic contact line ^(t) = max(h^^(£, x)) is indicated by the diamond shaped markers along the 
x-axis, equidistantly spaced in time. Note that the x-axis is compressed by a factor 1/e = 100. 


MATHEMATICAL MODEL 
Elasticity 

We begin with a description of the elasticity of a long thin and very wide sheet of length L [m], thickness b[m] ^ L, 
separated from a solid substrate by a viscous thin film of height h = h{x^t) [m] and h ^ where we denote the 
dimensions of the variables within [•]. The elastic sheet is attracted to a substrate by an intermolecular adhesion 

potential (Fig. 1 ), consisting of a van der Waals attraction and a strong short range-repulsion ^ 

[29] . with A [J] the Hamaker constant. Here a [m] is the equilibrium height where the attraction and the repulsion 
pressure balance; evaluating the energy at this minimum yields the effective solid-solid surface tension 7 = ^(cr) = 
A/ 8 cr^ [J/m^]. Then the shape of the sheet is determined by the balance of transverse forces on the sheet and yields 
the following equation for long wavelength elastic deformations m 

P = Bhxxxx ^h‘ ( 1 ) 

Here p is the local pressure along the sheet, B [J] is its bending stiffness and {')x = d{-)/dx^ and we have further 
assumed that <C 1 , so that we can neglect any geometrical nonlinearities in the shape of the sheet. 

For a static sheet, the pressure p = 0. In this case, multiplying 0 by hx and integrating by parts, or alternatively 
using the principle of virtual work, yields the boundary conditions at the static elastic contact line h{R) = hx{R) = 
0,hxx{B) = derived and experimentally verified nearly a century ago by Obreimoff [28], see Supplemental 

Information (SI). 


Hydrodynamics 

When the sheet is allowed to move in response to the attractive adhesive interaction which brings it into proximity 
with the substrate, the interstitial fluid begins to be squeezed out. On small scales, the inertia of the fluid and the 
sheet is dominated by viscous forces in the gap. In low Reynolds number limit, i.e. Re = ph‘^c/{Lfji) 1 with p 
[kg/m^] the liquid density, hoo the initial film height, c [m/s] the adhesion speed along the substrate, and p [Pa- s] the 
fluid viscosity, the fluid flow field, u = {u{x,yR),v{x,yR)) is governed by the incompressible two-dimensional Stokes 
flow equations Vp = /iAu, V • u = 0. Evaluating the different terms in the Stokes equations show that to leading 
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order in e = h^/L 1 the thin film fiow can be described by the lubrication equations [31] 

Px “1“ p^yy — d? Py — d? — d. (2) 

Using the no-slip boundary conditions on the two solid surfaces 

i4(x,0)=0, u{x,h) = 0, 

'z;(x,0) = 0, v{x,h) = ht, (3) 

to integrate (j^ yields u{x^y) = ^ {y‘^ — yh). On substituting this result into the continuity equation together with 
0 yields a sixth order nonlinear partial differential equation for the thin film height [32] 

) = 0. (4) 



Scaling, non-dimensional numbers and boundary conditions 


Scaling Q using the following definitions for the dimensionless variables; 


r h ^ X ^ ph'i 


t = t/r^ = 


tA 


12/ih^ 


( 5 ) 


give three non-dimensional numbers: B = BlrA^/AL^^ the ratio between elastic bending pressure and adhesion 
pressure; a = {(jjhoo), the ratio between the equilibrium adhesion height and the initial film height; and the aspect 
ratio e = h^/L. All our results in the text are given in both dimensional and scaled forms, but our figures show only 
the non-dimensional results scaled according to 

In the context of experiments on wafer bonding^ [TO] with /i = 2 x 10“^ Pa-s, hoo = 100 nm, L = 10cm, A = 10“^^ 
N-m [29], cr G [12 — 24] nm, B = 10“^ — 10 N-m, we find that e = 10“^ ^ 1, d G [0.06,0.12], B G [10“^ — 10“^). 
Furthermore, since the bonding speed c G 0.1 — 2 cm/s Eiaiin], the Reynolds number Re = ph‘^c/(pL) « 10 ^ ^ 1 
justifying the lubrication approximation Q. We note that since B 1 and multiplies the highest order derivative 
in Q we expect a boundary layer to form near the apparent contact line. Indeed, the size of this boundary layer 
Ic can be estimated by balancing the bending pressure Bhxxxx ^ ^/it arid the attractive van der Waals pressure 
AjZln?^ as h ^ cr so that Ic ~ a{3B/A)^ and is the length of the inner adhesion zone (see SI). Furthermore, we note 
that the slope near the apparent contact line scales as hx ~ cr/Z^ ~ {A/3B)^ 1, justifying our long wavelength 

approximation. 

In order to complete the problem formulation, we need to prescribe six boundary conditions at the two edges of 
the sheet. Assuming that the sheet is free of forces, torques and exposed to an ambient fiuid environment at constant 
pressure, this implies that at either end. 


hxx\x=0,L — hxxx\x=0,L — P\x=0,L — 0 - ( 6 ) 

Our mathematical model provides a general and compact transient formulation of the entire elastohydrodynamic 
adhesion phenomenon that includes the microscopic physics, extending and complementing previous steady-state 
analyses 


ANALYSIS 

To understand the behavior of the solution of the initial boundary value problem we solve equation @ 11 ) numerically 
by using a second-order finite difference method with a spatial discretization 5x/L = 1/400 — 1/1600 and use a Gear 
method [33] for the adaptive-time marching. In Fig. we show the results obtained by solving @ numerically, 
delineating three regimes, starting with the onset of adhesion at short times, an intermediate time regime as the ffuid 
filled bump is formed and a long time regime associated with the steady motion of the adhesive front R{t) and the 
fiuid bump. 
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FIG. 3: Waiting time for onset of adhesion. The scaled waiting time = C{3B/A)i = for the onset of adhesion 

as a function the the scaled bending stiffness B G [5 x 10“^ — 10“^), the scaled cut-off height a G [0.03 — 0.12] and the aspect 
ratio is e = 10“^, with C = 4. Square shaped markers represent S G [5 x 10“^ — 10“^), d = 0.12 and the star shaped markers 
(they are coincident) represent B = 10“®,d = [0.03,0.06,0.12], the dashed line represents the non-dimensional form of the 
scaling law in 0 - 



FIG. 4: A snapshot of the propagating adhesion zone. The fluid hlled bump ahead of the adhesion zone moves with a speed 
dR{t)/dt = c. As the adhesion zone grows, it displaces the viscous fluid into the bump making its horizontal extent i{t) 
and height H(t) to grow in time. The circle denotes the small inner region where viscous, elastic and adhesive forces are all 
comparable in magnitude. Note that the x-axis is compressed by a factor 1/e — 100. 
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a) b) 

FIG. 5: Evolution of the sheet profile onto a self-similar shape, a) A rescaled version of Fig. [fusing the form f{r]) = h{x, 
with f] = (x — ct)/i{t), where H{t) and i{t) are given in |T^ . The shapes of the sheet collapses onto a universal self-similar 
shape, here plotted at twelve different points in time starting from left-to-right for t E [48 — 438] with B = 5 x 10“^, a = 0.12 
and the speed c = 0.0133 is measured numerically, b) One of the shapes of the sheet (solid line) for t — 438 in logarithmic 
coordinates, which shows that the slope of f{inj) is well matched with the analytical prediction (dashed line) in the self-similar 
regime f(rf) ~ 773 . 


Waiting time 

As the attractive van der Waals pressure brings the surfaces together, fluid is squeezed out. The time for the 
onset of adhesive contact i.e. ^ a is associated with the time to displace the viscous fluid beneath the adhesion 
zone. Balancing all three terms in @ yields the horizontal scale over which pressure and adhesion balance each other 
when the height of the film is /loo, be. loo ~ hoo{B/A)^, and thence a scaling law for the viscous drainage time 

Tw (7) 

where (7 is a dimensionless constant. In dimensionless form this reads as . Our simulations 

of [4]-[6] follow this scaling prediction 0 . with (7^4 see Fig. It is useful to contrast this relation with the 
capillary (tension) waiting-time for thin fluid film rupture « ^Fere is the liquid-vapor surface 

tension, showing a weaker dependence on both A and hoo- 


Self-similar shape 

After first contact when h = a, the adhesion zone grows by displacing the fluid into a localized bump of height H{t) 
and width i{t)^ as shown in Fig. and appears to converge to a self-similar traveling form. We also see the formation 
of an apparent contact line, which we define as the position along the sheet where the curvature is a maximum i.e. 
R{t) = max(/ia,a,(x, t)). In the vicinity of the apparent contact line max(/ia,a,) viscous, adhesion and bending forces are 
all of the same magnitude (see SI). However, as we move away from the apparent contact line the adhesion pressure 
becomes small and the bump shape is described by elastic bending and fluid conservation, and constrained by the 
far field boundary conditions. To understand these features, we now turn to a simplified analytical description of the 
shape and speed of the adhesion zone and also determine the effective boundary conditions at the apparent contact 
line. 

We start by looking for the asymptotic form of the bump ahead of the adhesion zone, where adhesive effects may 
be neglected locally. We use a self-similar ansatz seeking a solution to 0 for the film height and the fluid pressure 
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FIG. 6: Speed of the apparent elastohydrodynamic contact line. The dimensionless adhesion front speed c as a function of the 
scaled initial condition corresponding to flat, linear and quadratic initial shape i.e. , t = 0) = [1 — 0.1 x (1 — tanh(25x)^, 0.5 + 
2x/3, l/2+2:r^/27], respectively. The scaled bending stiffness B G [5 x 10“^ —10“^] and the scaled cut-off height a G [0.03 — 0.12]. 
Squares - flat, d = 0.12, Be [5 x 10“^ — 10“^]; circles - linear d = 0.12, B G [5 x 10“^ — 3 x 10“^]; diamonds - quadratic, 
d = 0.12, B E [b X 10“^ — 7 X 10“®]; stars - flat, d G [0.03,0.06,0.12], B = 10“®. The dashed line represents the scaled form 
of IFl i-e. C = 3eV(40<7^Bi). 


of the form, 


h{x,t) = H{t)f{ri), 7] = 
p{x,t) = P{t)g{r]). 


X — ct 
£{t) 


(8) 


Here (•)' = d{')/dr]^ where r] is the scaled coordinate frame of the apparent contact line moving at the speed c, f{r]) 
is the scaled film height, g{r]) is the scaled pressure, and H{t), £{t) and P{t) are three unknown scaling parameters 
as illustrated in Fig. 

In the regime associated with the steadily moving adhesion zone, the characteristic height of the film H{t) ^ a so 
that the adhesive pressure A/3H^{t)f^{r]) <C 1 and the total pressure p(x, t) ^ Bh^xxx ~ Furthermore, 

ht = —cf'{r])H{t)/i{t)Pcf{r])dH{t)/dt—f'{r])H{t)r]di{t)/dt/i{t) ^ — c/'( 77 )^ 7 (t)/-^(t), where the first term is dominant 
at long times [35] . Then, on substituting the form into Q we obtain the asymptotic equation 

_ H\t)Bp{g)r'{g) ^ 

12gcP{t) ■ ^ ^ 

Defining the pre-factor k = H^{t)B/12 (|^ reads as —kp{r])f""'{r]) = 1 with a polynomial solution given by 


fid) = g{v) = 


37 ] 


k = 


BH^{t) 


243 


( 10 ) 


7k ’ 12]icP{t) 

Having obtained the form of the solution /(r^), we now turn to determine the rate of fluid mass swept out by the 
bump and thence determine its vertical and horizontal scale factors H{t) and i{t). For a film with an initial constant 
height d(x, t = 0) = the volume swept by the bump is given by 


net pOO 

V= h^odx = hooct = H{t)£{t) / f{ri)dri. 
J 0 J — 00 


( 11 ) 
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We assume that the scaled displaced mass is normalized using the condition f{ri)dri « f{ri)dri = 1, consistent 
with the numerically measured mass in the bump fa 0.85. Then, from (|10|) and^|ll|) we find that 


H{t) 




12iik ) 



( 12 ) 


i.e. the shape of the fluid bump ahead of the adhesion zone is a function of the properties of the elastic sheet, the 
fluid, and the speed of the adhesion front. In Fig. we corroborate these scaling predictions by showing that our 
numerical simulations of Q collapse to a universal rescaled shape. Our numerical observations of the bump ahead of 
the adhesion zone are consistent with experimental observations m- Moreover our similarity solution reveals that in 
the immediate vicinity of the steadily moving front h{x^t) ^ ^ H{t) {{x — ct)/£{t))^ (Fig. |^), consistent 

with experimental observations near the adhesion zone [ 8 ] . While we have focused here on the case when the far field 
corresponds to a flat sheet, we can generalize our analysis to an arbitrary initial condition h( 0 ,x) oc x^, so that the 
volume swept up by the bump depends on n via the shape of the sheet in the far field (see SI). 


Adhesion speed 


Having characterized the self-similar shape of the traveling bump, we now determine the speed c of the apparent 
contact line. In a traveling wave frame moving with a speed c, our numerical simulations show that the elastic sheet 
squeezes the liquid ahead, forming a bump that grows with time as it sweeps up the fluid ahead. Once the contact 
zone reaches a steady state, we expect that no energy is expended in local bending of the sheet, since the size of 
the contact zone Ic is invariant. Thus, the speed of the apparent contact line should be determined by the balance 
between the driving adhesive power and the viscous dissipation in the contact zone. 

At a scaling level, the adhesive power per unit width H^ scales as the product of the speed and the energy per unit 
area of adhesion, i.e. H^ ^ c^(cr) ^ cy ^ cA/a^, while the viscous power per unit width H^^ scales as /i(Vt/)^f^, where 
VI/ ^ c/cr, and the volume per unit width ft ^ alc^ so that H^^ ^ fic^ldcf. Balancing the adhesive pover with the 
dissipation i.e. H^ ^ H^^ yields c ^ A 4 ^ which we may rewrite as c ^ 7 A 4 ^ an expression very similar 

to that given in [ 8 ], although we note that by accounting for microscopic adhesion, we do not need to postulate the 
existence of an apparent slip length used in [ 8 ]. 

Going beyond this simple scaling analysis, we write the exact balance between the dissipation rate and the adhesion 
and bending power in the steady state problem as [ 8 ] 


nOO 

Jo 


—^dx = — 
12 /i dt 


f 


+ ) dx 


(13) 


To understand the relative magnitudes of the different terms, we 


introduce the similarity variables ( 8|12 ) into (13). 


In the contact zone, the elastic bending pressure and the attractive van der Waals pressure balance each other so that 
^ = Bhxxxx^ he. ~ approach the apparent contact line, where the sheet is statically 

adherent, with r] ^ 6{t) and / = 77 ?, we find that 


~ /"TV = f ^ 

^ ^ \3B) \H{t)) \36iic'^kho 

First, we note that the rate of change of elastic bending energy reads as 

1 /•! 

{B^{l2^ikfhooC^Y t-i / r]-idr]. 

Js(t) 


_ 3 

t 8 . 


(14) 


(15) 


Inserting (12|l4) in (15) shows that at its lower limit the integral scales as t 8 ( 5 (t )3 t ^ and at the upper limit 

1__l^ 

the integral scales as t“ 8 . Thus, although the integrand diverges in the vicinity of the contact line, the integral itself 
becomes vanishingly small at long times. Physically this is consistent with the idea that bending energy is neither 
stored nor dissipated as the front has an invariant shape that persists with time. 
















FIG. 7: Characterizing the elastohydrodynamic contact zone. The scaled maximum curvature in the contact zone max(B 2 hxx) 
with the scaled bending stiffness B — and three different scaled minimum heights a — [0.03,0.06,0.12]. The dashed line 
is the scaled form of the relation (19) and illustrates the singular nature of the curvature as d ^ 0. Inset: The curvature of 
the sheet hxx with a = 0.06, where the star shaped marker dehnes the position where the curvature is a maximum. 


Next, we note that the viscous dissipation rate reads as 


f 


h^Pl 

12/u 


dx ^ — 


Jim 


f 

Js(t] 


T] 3 dr]. 


(16) 


In the neighborhood of the contact line, the integral remains bounded due to the presence of the cut-off length 6{t). 
Indeed, in the asymptotic regime t ^ oc the value of the integral at its upper limit scales as oc and thus becomes 


vanishingly small. The value of the integral is instead dominated by the lower limit. Inserting 6{t) from (14) in (16) 
yields an estimate for the viscous dissipation as ^ jA)^ ^ which is constant for a front traveling with a constant 

speed, consistent with our simulations (see SI). 

Finally, we note that the adhesive power is dominated by the contribution from the maximum of the adhesion 
potential i.e. 


f] poo j pR{t) 

— / ^{h)dx ^ — / max (ff>(d)) dx 
dt Jq dt Jq 


Ac 


(17) 


Balancing the rate of change of adhesion energy 0 with the viscous dissipation rate ( p^ yields a scaling law for the 
adhesion front speed 


A 5 

^44 

160 a‘^pLBi ’ 


(18) 


which in dimensionless form reads as c = 3e^/(40(j^55), consistent with our earlier simple scaling estimate. Since 
A/8cr^ = 7 , the interfacial tension, we see that the speed is independent of the microscopic length scale cr, but does 
depend on the adhesion strength A, consistent with earlier results [8]. 

In Fig. we show the speed c of the adhesion front obtained by solving for a range of parameter values of 
B e [5 X 10 ^ — 10 ^], initial far-field shapes determined by n, cut-off heights a G [0.03 — 0.12] for the aspect ratio 
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e = 10“^, and see that the simulations corroborate (18). Our analytical prediction also favorably compares with 
experiments (SUHUTO]. By using their material parameters in (18) i. e. 7 = 0.12 N/m, /r = 2 x 10“^ Pa-s, B = A 
Nm [HI [To] and A = 3 x 10 “^^ Nm [29], gives a theoretical wafer bonding speed c ^ 0.5 cm/s consistent with the 
experimental observations that yield c ~ 0.2 — 2cm/g[5l l8l ITO]. 


Boundary conditions at the apparent elastohydrodynamic contact line 

We now turn to addres s the question of the singular nature of the shape of the adherent sheet in the adhesion zone 
and by using ([^ [l^|l2|l4[ ) we can evaluate the similarity form for hx ^ ^ H^ H(t)5i /l{t) « 

(^)4. Thus, our analysis suggests that the apparent contact angle hx of the elastic film approaches a constant value 
in the contact zone, similar to the Young-Laplace contact condition for a static fluid-fluid contact line on a solid 
substrate [25H27] . This implies that the curvature of the sheet diverges in this zone. Of course, this divergence is 
suppressed in the adhesion zone over a size Ic ~ g{B jA)^ determined by the balance between the bending and the 
attractive van der Waals pressure. Indeed, in the contact zone hx ~ crjlc ~ ^xx ~ cr/^c ^ {A/. 

In Fig. we show that the scaled maximum curvature in the contact zone max(ha,a,) is in agreement with both this 
scaling law and the similarity solution. 

Finally, we determine how flow and elastic deformation processes together yield the effective boundary conditions 
associated with mesoscopic deformations of the elastic sheet, viewed from the perspective of the outer problem where 
the microscopic adhesion is relatively unimportant. Approaching the dynamic apparent elastic contact line R(t) from 
the right side, we note that outside the contact zone, the van der Waals adhesion pressure is negligible, so that the 
equation of motion simplifies to 12/iht — {Bh?hxxxxx)x = 0- Looking for a traveling wave solution of the form 
h{xA) = g{x — ct) then simplifies the equation to the form —eg' — {Bg^g'"")' = 0, where (•)' = d{')/d{x — ct) which 
has a similarity solution of the form g ^ [x — consistent with [8], and our earlier similarity analysis. We note 

that as we approach the contact line, x — ct ^ 0, ^ ^ 0, ^ 0, but the curvature g" oc. However, this divergence 

is ameliorated by the presence of the inner scale associated with the contact zone where microscopic adhesion, bending 
and fluid flow are in balance. This leads us to the conclusion that the effective boundary conditions for the outer 
problem at the contact line i?(t), i.e. ignoring the details of the contact zone, are given as 

/i(i?)=0, hJR) = 0, c=^=Ki -——r (19) 

where Ki = 0.8 ±0.1 is a dimensionless constant evaluated numerically. In dimensionless form these read as; 
h{R) = 0Ax{R) = 0, c = /{AOa^B^). It is useful to contrast these conditions with those for a static adhered 

sheet, where at the elastic contact line [28] i.e. h{R) = hx{R) = 0, hxx{R) = static condition 

for the curvature at the contact line is replaced by a condition for the speed of the apparent contact line. In the limit 
cr = 0, A = 0, the boundary conditions in ( [l^ read h = hx = c = 0, and the contact line becomes stationary da El], 
consistent with the classical boundary conditions for a non-adherent elastic sheet in contact with a solid. 

These results might also be contrasted with the analogous problem of describing static and dynamic contact lines 
in interfacial hydrodynamics where a static contact line has a constant contact angle, whereas a dynamic contact line 
has a contact angle that is a function of its speed [25H27] . 


CONCLUSIONS 

Our theoretical study of viscously limited elastohydrodynamics of adherent sheets captures the entire process of 
dynamical adhesion, from the short-time onset of adhesion to the long time dynamics associated with a steadily 
propagating adhesion front. A simple mathematical model provides a compact formulation that naturally couples 
the microscopic physics at the apparent elastic contact line and the macroscopic physics associated with elastic 
deformation and fluid flow, complementing earlier preliminary analyses of the problem. Numerical simulations reveal 
different regions of the sheet; an adherent zone x < R{t) at constant height a, an inner contact zone where viscous flow, 
elastic bending and microscale adhesion are all equally important, and an outer region with h^ a where the sheet is 
described by viscous flow and elastic bending. An asymptotic similarity analysis of the governing partial differential 
equation allows us to describe the outer zone consistent with the solutions of the governing partial differential equation, 
and leads to a self-similar shape of the elastic sheet and its propagation speed. In addition, we have derived the 
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effective boundary conditions for the dynamic apparent elastic contact line, which highlights its singular nature and 
distinguishes it from analogous contact line conditions for a static sheet. 

Just as high resolution Total Internal Reflection Fluorescence microscopy allowed for a more detailed view of the 
dynamics of liquid contact lines [27] , we hope that our theoretical study might engender further investigation of the 
nature of the elastohydrodynamic contact line that arises in a wide range of problems where elastic interface adheres 
to solid substrates in fluid environments. 
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Fig. A 1: Schematic. A static elastic sheet that is adhered to a solid substrate at is left edge at an equilibrium height h = cr, 
while the right end is supported at a height hoo- Note that the x-axis is compressed by a factor 1/e = 100. 


APPENDIX 
Static sheet 

Long wavelength deflections in a static sheet yield the equation of equilibrium for its shape that is given by 

P = Bhxxxx 

The accompa nying boundary conditions are /i(0) = hx{0) = 0 and at the support, h{L) = hoo, hxx{L) = 0, see Fig. % 
We note that Jai| follows as the Euler-Lagrange equation associated with the energy functional f^[Bh‘l^/2-\-^{h)]dx. 
Taking the variation of this energy with respect to the location of the the contact line at x = R, the principle of 
virtual work yields the contact line conditions: 

h{R) = cr, hx{R) = 0, Bhl^{R) = 2^(cr), (A2) 

where the last condition can be rewritten as hxx{R) = ( 4 ^^)^ ~ Similarly, these boundary conditions can be 

obtained by multiplying [1] by hx and integrating by parts. 


Role of initial shape of the sheet on its dynamics 

Different initial shapes affect the dynamic shape of the sheet during adhesion. First we recall the similarity ansatz 
and the polynomial solution to the self-similar equation in the regime where A/h^ ^ 1, i.e. where adhesion can be 
neglected. 


h{x,t) = H{t)f{r]), T] = 

fiv) aiv) = 

We account for a generic initial shape of the form h{x, t = 0) « ehoo{xlL)'^ corresponding to flat, linear and quadratic 
profiles and e is a constant pre-factor. The volume swept by the bump is given by 

V(n) « r eh^i^rdx « ^^ooe(cyL)»+^ ^ ^ ^ 

Jo L n + 1 J_^ 

We normalize the scaled displaced mass using the condition f{r])dr] ~ f{r])dr] = 1 and the numerical simu¬ 
lations indicate that this is a fairly good approximation, as seen in Fig. 5 where the numerical area of the bump is 
0.85. 


X — ct 

TT’ 


k = 


_ 243 


(A3) 
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a) 


b) 


Fig. A 2: Self-similar shape of the adhesion front, a) The scaled shape of the sheet h for three different initial conditions 
n = [0,1,2] for t = [438,458,478], respectively. h{x,t = 0) is; constant height n = 0, h{x,t = 0) = 1 — 0.1 x (1 — tanh(25x)^, 
linear height n = 1, h{x, t = 0 ) = 1/2 + 2x/3, quadratic height n = 2, h{x,t = 0) = 1/2 + 2f^/27. b) Rescaled height and 
x-coordinate according to ( |A5| ), i.e. f{r]) = r) — {x — ct)ll{t) gives the self-similar shape of the bump.We need different 

pre-factors to collapse the shapes i.e n = 1, 0.8 x [H(t^n — l),^(t,n = 1)] andn = 2, 0.66 x [H(t^n — 2),^(t,n = 2)], however 
both pre-factors are 0 ( 1 ) and the influence on the detailed mass balance / ( 77)^77 may come from the far-fie Id shape of the 
sheet that requires matching to the self-similar region. 


Then, substituting ( |A3| ) in ( |A4[ ) we find that 

( 




r BV^{n) Y 

( 12 /Ltcfc ) ’ 


(A5) 


i.e. the dynamic shape of the elastic sheet is a function of the fluid, the adhesion speed c and the initial film shape 
through V{n). We see that the initial condition determines the dynamics of the height and width of the bump, 
corresponding to; constant height [n = 0, H{t) oc ts, oc ts], linear height, [n = 1, H{t) oc t 4 , l{t) oc t^] and 
quadratic height [n = 2, H{t) oc ts, oc t“^] but does not affect the velocity of the adhesion front, see Fig. 6 . To 
test if the self-similar ansatz holds we have scaled three membrane shapes at different points in time and with very 
different mass accumulated in the bump, see Fig. Although the dynamics of the bump is a function of the initial 

condition, it appears to converge to a universal self-similar shape. Note that a pre-factor [0.8, 0.66] for both H{t) and 
i{t) is adjusted for the scaling of n = [1,2], respectively. Determining this scaling factor analytically requires us to 
match the different regions which is a function of the initial condition that appears to have a small influence on the 
volume integral / ( 77 )^ 77 . 


Dissipation 

Since the apparent elastic contact line is found to propagate with a constant velocity, we anticipate that the viscous 
dissipation rate in the thin fluid film is constant. To verify this, we extract the viscous dissipation rate, given by 
equation [16] in the main text, from one of our numerical simulations as shown in Fig. % it approaches a constant 
value after the initial transients die out. 
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Fig. A 3: The solid line shows the non-dimensional viscous dissipation as a function of time for the parameters; B = 3 x 10~®, 
n = 0, e(0) = 1 d = 0.12. At short times, as the adhesion zone is formed, viscous dissipation rapidly increases. At long times 
the adhesion front propagates with a constant speed c and the viscous dissipation approaches a constant value. The dashed 
line is intended as a guide to the eye. 






